Long-lived memory for electronic spin in a quantum dot: Numerical analysis 
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Techniques for coherent control of electron spin-nuclear spin interactions in quantum dots can 
be directly applied in spintronics and in quantum information processing. In this work we study 
numerically the interaction of electron and nuclear spins in the context of storing the spin-state 
of an electron in a collective state of nuclear spins. We take into account the errors inherent 
in a realistic system: the incomplete polarization of the bath of nuclear spins and the different 
hyperfine interactions between the electron and individual nuclei in the quantum dot. Although 
these imperfections deteriorate the fidelity of the quantum information retrieval, we find reasonable 
fidelities are achievable for modest bath polarizations. 
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I. INTRODUCTION 



Currently, the development of techniques for coher- 
ent manipulation and control of electron spins in semi- 
conductor nanostructures is an area of vibrant re- 
search activitji 2 * 3 ^*^. Such techniques are needed for 
many coherent spintronic devices^, such as a Datta-Das 
transistor—, and for quantum computing using the spin of 
an electron in a quantum dot (QD) as a qubiti^>. How- 
ever, the decoherence time of an electron spin constitutes 
a fundamental obstacle on the way to realization of a QD- 
bascd quantum computation 9,10 . In order to enhance the 
coherence time, a quantum memory protocol has been re- 
cently suggested 11 which allows mapping of the quantum 
state of the electron spin on the collective spin of the nu- 
clei present in the QD, thereby taking advantage of the 
intrinsically longer coherence times of nuclear spins. Un- 
der ideal circumstances, the proposed scheme allows to 
encode and, at later time, to retrieve, the electron spin 
state with near 100% reliability. Furthermore, it can 
be extended to implementing arbitrary 1- and 2-qubit 
operations^. Imperfections, which should be expected 
in realistic settings, lead to errors, thus deteriorating the 
fidelity of the quantum memory protocol 11,13 . Some er- 
rors can be reducecU^i^, but some are inherent to the 
protocol, such as the incomplete polarization of the nu- 
clear spins and the spread in the hyperfine interactions 
between the electron and the nuclei in the QD. In this 
work, we use numerical simulationsi£ii£ to analyze the 
effect of these errors on the performance of the quantum 
memory protocol. We show that the scheme remains fea- 
sible for the spin bath polarizations exceeding 80%, with 
the minimum fidelity of ~ 80%. 

First, let us recall the scheme of Ref. [H} The electron 
with a spinor wave function \<f>o) = a\ I) + j3\ |) is in- 
jected into an empty quantum dot, and starts interacting 
with the bath of N nuclear spins in the dot. The most 
important part of this interaction is the isotropic con- 
tact hyperfine coupling. The total Hamiltonian of the 
system, taking into account the Zeeman energy of the 



electron spin in the external field Hq, becomes 

N 

H = g^BH S z +Y / ^klkS (1) 

fc=i 

where 5=1/2 is the electron spin, Ik — 3/2 are the spins 
of the nuclei (the fact that the nuclear spins of Ga and 
As are not 1/2 is not important here, up to renormal- 
ization of the appropriate parameters, and we will con- 
sider them as spin- 1/2 in the following discussion), and 
Ak = (87r/3)g e ^Bg„^„|\E'(xfc)| 2 is the contact coupling 
determined by the electron density |\&(xfc)| 2 at the site 
Xfe of A;-th bath spin, and by the Lande factors g e , g n and 
the magnetons [Ab, fJ-n of the electron and the nuclei, 
respectively. At maximum polarization, the nuclei can 
produce an effective field J^k Aklk/gl^B = 5.2 Teslai&. 

Other terms in the Hamiltonian can be omitted 
when discussing the storage and retrieval processes due 
to their smallness at the relevant timescales. These are, 
e.g., the Zeeman energy of the nuclear spins, the dipolar 
coupling between the bath spins, or between the electron 
spin and the bath spins. Also, we assume that the injec- 
tion is performed quickly, so that the approximation of 
the sudden change of the Hamiltonian is applicable, and 
neither the state of the bath nor the state of the electron 
spin are affected by the injection process. Alternatively, 
fast control of the Zeeman splitting of the electron allows 
for the electron spin-nuclear spin system to be brought 
into resonance in a sudden approximation, and is similar 
in nature to the controlled injection/ejection of electrons 
from the quantum dot. 

If the electronic wave function were spread uniformly 
over the dot volume, making all Ak equal to A, and if 
the nuclear spins were completely polarized (i.e., if the 
bath were in the state \xo) = \ I, l> ■ ■ ■ D) then, for the 
external field adjusted to be Ho — A(N — l)/(2g e fig), 
the motion of the compound system (the electron and 
the bath spins) would become a sinusoidal oscillation be- 
tween the states \4> ) <8> |xo) and | J.) ® [a|xo) + */?IXi)] 
where |xi) = (l/VN) J^fe S£\Xo) (N is the number of 
the nuclei in the dot). Thus, after half-period of such 
oscillations, the initial state of the electron spin becomes 
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encoded in the state of the bath, and the electron can 
be ejected from the quantum dot. In order to retrieve 
this information after some time, another electron in the 
state | J.) can be injected, and after a half-period of 
the oscillations, the state of the electron spin is again 
\4>q) = a\ I) + (3\ T)- Of course, for this two-step scheme 
to work, the delay between the information encoding and 
retrieval should be small in comparison with the charac- 
teristic times of the omitted terms in the Hamiltonian 
(in practice, these terms should be reduced using, e.g., 
NMR techniques*^). 



However, real baths are incompletely polarized and the 
approximation of complete polarization is never applica- 
ble to a real bath with N ~ 10 4 - 10 6 . To explain this, 
let us consider a highly polarized state of the bath de- 
scribed by a density matrix p = (1/Z) exp (— 7P), where 
P = ^2 k Ik, and Z is chosen such that Tip = 1. For 
such a bath, the polarization is P = tanli7/2, and for 
large 7 the value of AP = 1 — P is very small. How- 
ever, the statistical weight w(M) of the states with a 
given value M of P is w(M) = fe (l - 6) N ~ k , where 
9 = e7(l + ei), the integer k=N/2- M, and Cf is a 
binomial coefficient. The distribution w(M) is approxi- 
mately Gaussian centered at M — (—N/2)(l — AP) with 
the rms a 2 = (N/i)AP(2-AP). Therefore, unless AP is 
unrealistically small (AP <C 1/VW ~ 10~ 2 ), a very large 
number of the bath states are involved in the time evolu- 
tion, as opposed to only two states |xo) and |xi) involved 
in the case of Ak — A and AP = 0. All these states 
have different oscillations frequencies, and form complex 
time-dependent superpositions, so that the fidelity of the 
information encoding and retrieval is reduced. Therefore, 
the errors associated with AP > (incomplete polariza- 
tion) and Ak ^ A (spread in the hyperfine couplings) are 
intrinsic to the protocol. Curiously, the symmetry prop- 
erties of the spin interaction allows for some correction, 
by state engineering and spin-echo techniques, for these 

14 

errors**. 



The paper is organized as follows. In section II, we 
consider qualitatively the analytical behavior of the stor- 
age protocol when inhomogeneity is negligible (Ak — A) 
but polarization is less than perfect, and indicate how 
many of the effects of different oscillation frequencies and 
Hilbert space dimension lead to only moderate decreases 
in the fidelity of storage and retrieval. We then evaluate 
the same situation numerical in section III, and develop 
the useful measures of performance for the protocol. In 
section IV, the numerical results are extended to the case 
of inhomogeneous interaction (Ak 7^ Ak'), and then in 
section V the minimal fidelity of the protocol (minimized 
over all possible electron spin states) is calculated. 



ii. qualitative analytics for a h = a: 
why the protocol works for 

i/Vn< ap< 1 

As indicated abovei*, for Ak = A and AP > the 
fidelity of the quantum memory protocol remains high 
if AP < 1, even for AP > 1/VN. It is instructive to 
understand this robustness in detail. 

To grasp the idea, let us consider only the first, encod- 
ing, step. We note that for Ak = A, the total spin of 
the nuclear bath / is an integral of motion. If the initial 
state of the bath is \Iq,Mo), where io is the value of / 
and Mq is the value of P , then the time-dependent wave 
function of the compound system is 

m)) = /3|Vi(*)> + a|V> 2 (*)) (2) 

\rpi(t)) = (coswit-i£isinu;it)| T}®|/o,M ) 

- iiji sina>ii| |) (g) |I , M + 1) 
\tp 2 (t)) = (cosw2* + i6sinco 2 *)| 1)®\I Q ,M Q ) 

- ir] 2 sinu 2 t\ f) <E> \Io,M - 1) 

where io± — y/hf + A 2 Cf/4, £1 = hi/iJi, and 771 = 
AC\/(2loi) (lu 2 , £2 and i] 2 are defined similarly). The 
relevant matrix elements are 

h x = g e p, B H /2 + A(M + l/2)/2, 

h 2 = g e ^ B H /2 + A(M - l/2)/2 , 

Ci = ^(/o-AfoXJo + Mo + l) , 

C 2 = V(/ +Af )(Jo-Mo + l) ■ 

The oscillation frequencies u>i j2 are strongly dependent 
on both Jo and Mq. For the bath's initial density matrix 
p = (1/Z) exp (— 7P), the statistical weight for the state 
|ib, Mo) is 

w(I , M ) = w(M Q )(C% - C^/Cf (3) 

where k — N/2 — Mo and m = N/2 — Iq are integers, 
and is a binomial coefficient. For AP <C 1, both 
N — k and m are small (note that Mq < 0, so that k is 
close to N), and for a given Mq, it is easy to see that 
w(I ,M ) ~ AP {N - k ~ m \ The main contribution to the 
evolution comes only from the states with Io = —Mq, and 
the states with larger values of Io(> —Mo) have too small 
a weight to impact the oscillations, and, correspondingly, 
to reduce the fidelity. Thus, the distribution of the fre- 
quencies oj\ i2 in Eq.|5]is determined only by w(Mq), and a 
spread of the oscillation frequencies a ~ V NAP is much 
smaller than the mean frequency Q ~ v^V- 

Using Eqs.EEJit is easy to obtain an analytical approx- 
imation for N » 1 and 1/yN <C AP <C 1 in the leading 
order in AP (but uniformly in time), since £,i,rj 2 ~ 0, 
£2,771 ~ 1- For the spin initially polarized along the x- 
axis (a = (i = l/\/2), 

2 2 

s x (t) = cos ^N(\- AP)) exp + O(AP) 

(4) 
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while for the initial "spin- up" state, 

s z (t) = cos (^y/N(l- AP)) + O(AP) (5) 

where r = At, and only s x (t) is influenced in the lead- 
ing order in AP. The difference between the results for 
transverse (x) and parallel (z) states indicates that the 
memory protocol is not isotropic: its performance de- 
pends in part upon the spin state being stored. As such, 
in analyzing errors, we will consider several different per- 
formance measures, described in detail below. We first 
consider the reliability of numerical simulations. 

III. QUANTITATIVE ANALYTICS FOR A k = A 
AND RELIABILITY OF COMPUTATIONS FOR 
SMALL TV 

We focus on the case of Ak — A in order to assess 
reliability of our simulations with N = 20 (presented 
below, performed for the case of Ak ^ A) , and to quantify 
the fidelity in detail, without omission of higher-order 
terms in AP. The case of Ak — A is a good check since it 
allows analytical solution, so that the results for N = 20 
can be directly compared with the results for realistic 
value N ~ 10 4 . 

Let us focus first on the encoding step to determine 
what will be useful performance measures for the proto- 
col. Since hi ^ hi in Eqs.|2| the components S x , S y of 
the electron spin rotate with time in the x-y plane with 
the angular velocity of order of A/4. This rotation is 
uniform, and does not depend on the particular values 
s x , s y of S x and S v , because the Hamiltonian JIJ is in- 
variant with respect to the rotation in the x-y plane {Ti. 
commutes with exp [— ie(S z + I z )} for any e). As a result, 
even in perfect situation of AP = 0, A^ = A, for the ini- 
tial value of Sy = 0, the final retrieved value Sy ^ 0, 
and Sx < Sx . However, this error is noticeable only for 
relatively small N (N < 100), since the total time of the 
interaction of the electron spin with the bath is only one 
oscillation period, i.e., roughly 1/(A\/~N), while the an- 
gular velocity is of order of A/ A. Furthermore, this error 
is easily recoverable by application of a correcting pulse 
of the external field H z . Therefore, a physically mean- 
ingful measure of errors associated with decoherence of 
the S x operator is the total transverse component of the 
electron spin, 

s T = yj sl + si , (6) 

which is proportional to s£ if Sy = 0. This is illus- 
trated in Figs.^a-d), where we plotted the oscillations 
of s x and st for AP = (left column, (a) and (c)), and 
for AP = 0.2 (right column, (b) and (d)), for the initial 
condition Sy = s z — 0, si 1 ' = 1. The solid lines cor- 
respond to N — 20, and the dashed lines correspond to 
N = 10 4 . It is important that during the first oscillation 




(g) Time (h) 



FIG. 1: First step of the protocol, encoding. Oscillations are 
shown for s x (a,b), st (c,d), s z (e,f) and so (g,h), obtained 
for AP = (a,c,e,g) and for AP = 0.2 (b,d,f,h), with N = 20 
(solid curves) and N = 10 4 (dashed curves) . 



period, the results for st obtained with N = 20 are close 
to those obtained with N = 10 4 , so that N = 20 is suf- 
ficient to reliably estimate (with the precision of ~ 0.05) 
the fidelity of the quantum memory protocol at much 
larger values oi N. st is one of the three performance 
measures we use for the remainder of the paper. 

In addition to the transverse errors, the oscillations of 

the value s z of S z are not simply proportional to the 

initial value s^ , so that the retrieved final value si^ — 

Ui + U2S^ (i.e., the value s z ^ is linearly dependent on 

Sz but not proportional to it). This type of error is, as 

far as we know, irrecoverable. To quantify it, we consider 

the value of si for the case s z — 1, s£ = Sy = 0, and 

the value of ag for the initial conditions Sx = 1, Sy — 

s z 1 ' — 0. For the latter case, the value s z ^ is denoted 

(f) 

as so, to distinguish it from the value s z obtained with 
the initial conditions s^ — 1, s x ^ — s y ^ = 0. s z and 
sq are the remaining two performance measures used in 
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FIG. 2: Final step of the protocol, retrieval. The values so 
(solid line), s z (dashed line) and St (dotted line) are shown, 
obtained for AP = 0.2 with N = 20 (a) and N = 10 4 (b). 

the paper. In ideal situation of AP = 0, Ak — A, the 
value so should be zero, while in a non-ideal situation, sq 
quantifies the protocol error. Figs, me-h) illustrate this 
point, showing the oscillations of s z and sq for AP = 
(left column, (e) and (g)), and for AP = 0.2 (right 
column, (f) and (h)). The solid lines correspond to N = 
20, and the dashed lines correspond to N — 10 4 . Again, 
during the first oscillation period, the results obtained 
with N — 20 are close to those obtained with N = 10 4 , so 
that N = 20 is sufficient to reliably estimate the fidelity 
of the quantum memory protocol. 

Now we study performance of the quantum memory 
protocol at a quantitative level. The goal is to ensure 
that consideration of N = 20 is sufficient for making 
conclusions about the realistic case of N ~ 10 4 -10 6 . For 
Ak = A, using the Eqs. I2l.il the evolution of the reduced 
density matrix of the electron spin can be computed for 
the whole two-step protocol for N < 10 8 . The spin ejec- 
tion is equivalent to a measurement of S z described as a 
von Neumann's projection, so that at the ejection time 
t = t e , the wave function 

m e )) = | T) ® h) + 11)® h) (7) 

is transformed into a mixed state with the compound 
system's density matrix 

Re = I TXT I ® h)M + 1 0(1 1 ® h)H . (8) 

At the retrieval step, after injection of another electron 
in a dot, the compound system's density matrix becomes 

^ = u)ai®(ki)^ii + h)(^i) . (9) 

The system's subsequent evolution can be again calcu- 
lated using Eqs. |21 The resulting analytical expressions 
are complicated, but they can be easily estimated numer- 
ically for N < 10 s . 

Fig. |3 illustrates the final step of the protocol, the 
retrieval of the initial state of the electronic spin, for 
AP = 0.2 for N=20 (a) and N = 10 4 (b). All three pa- 
rameters quantifying the performance, st, s z and sq, are 
shown for times close to the time of the state retrieval 
(half-period after the injection of the second electron), 
and one can see that the curves for N — 20 almost ex- 
actly replicate the curves for N — 10 4 (the difference is 



~ 0.05). Moreover, it is important to note that the max- 
imum performance of the retrieval is achieved at slightly 
different times for sy, s z and sq. Below, we assume that 
the information retrieval is performed when s z reaches 
the maximum. In section llVl we evaluate all three per- 
formance measures, s z ,st, and So, while calculation of 
the fidelity of the protocol for arbitrary quantum states 
is considered in section Ivl 

IV. NUMERICS FOR A k / A: COMPUTATIONS 
FOR N = 20 AND ANALYSIS 

In a real quantum dot, the electron density is non- 
uniform, so that the values of the hyperfine coupling 
constants Ak differ substantially across the dot. For- 
mally, the Hamiltonian Q can be solved analytically for 
Ak 7^ A by Bethe ansatz 21 . However, this formal solution 
is impractical, especially in our problem, where the state 
of the bath has a complex form, with a large number of 
excited states, and where the non-unitary projection (see 
Eqs. 18191) is an important part of the protocol. We are 
not aware of any way of employing the formally exact 
solution of Ref. |2l] (which requires solving a system of 
order of 2 N /yN nonlinear equations) for practical quan- 
titative estimates in such complex problems as the one 
considered here, even for a small number of the bath 
spins. For similar reasons, numerical approaches based 
on the P-representation of the spin density matriji 2 ^ are 
also not applicable. Instead, we use the numerical ap- 
proach suggested in Ref. ITaTlTt the direct solution of the 
time-dependent Schrodinger equation for the compound 
system (the electron spin plus the bath spins). This al- 
lows straightforward numerical modeling of the quantum 
memory protocol. The drawback of such an approach is 
that the exponential scaling of computation costs with 
N limits simulations to N ~ 20 spins. Considering the 
above case of Ak = A, we have shown that N = 20 is suf- 
ficient to make our simulations reflect larger N behavior 
with a precision of order 0.05. 

In order to understand whether this remains true for 
Ak =/= const, let us consider the case of perfect polariza- 
tion P = 1 which can be studied analytically with preci- 
sion of 1/N for large N. First, let us consider the electron 
initially in the state \(f>o) = | T) and the bath in the state 
Ixo) = | Jo -L) • • •)■ By acting with the Hamiltonian on 
the state \ipo) ® |Xo)j we nn d that 

H\^o) = (g e » B /2)(H + H Ovh )\Tp ) + (&/2)|V>i) (10) 

where Hovh — ~^2k Ak/(2gelJ-B) is the Overhauser field 
acting from the nuclei on the electron spin, b — 
and = ||)® |xi) = 11)® [(1/VAOE^Ixo)]. 
Analogous result for \ipi) can be presented in the form 

HWl = -(ffeW2)(ffo + Ho«fc + M3/Ma)|Vi}(ll) 
+ (&/2)|Vo> + (&/2)K) 

where M n = ^ fc A^ is the n-th moment of the distri- 
bution of the hyperfine couplings Ak (correspondingly, 
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b = y/M2)i and the extra state \u') is orthogonal to both 
|^o ) and By choosing the external field equal to its 
optimal value 



-H 0vh - M 3 /(2M 2 g etJ , B ) 



(12) 



and subtracting a constant from the Hamiltonian, we 
have 



(6/2)1-0!) 
(6/2)|Vo> 



(13) 



(b/2)\u') 



where H = H + M 3 /(2M 2 ). 

This system is "practically" closed, since the norm 
of the extra state \v!) is small for large TV, \\u'\\ 2 = 
[(M 4 /M|) - (Ml /Ml)} ~ 1/TV for a "non-pathological" 
distribution of Afe. As a result, the calculations for both 
the encoding and the retrieval steps can be done explic- 
itly with the precision of l/N 2 . Introducing new basis 

states |</> + ) = (l/V2)[|Vi> + |Vo>], |^-> = (l/V2)[|Vi> - 
l^o)]) and |u'/||it'||), the Hamiltonian Ti can be repre- 
sented by the matrix 



1 5 
H = (b/2) | -1 -S 
S -5 



(14) 



where S = \\u'\\/y/2. The unitary evolution during the 
encoding and the retrieval steps is described by the evolu- 
tion operator exp (—iHt) which can be easily calculated, 
as well as the system's evolution during the projection 
step. We also note that the optimal times for encod- 
ing and for retrieval differ slightly (by the terms of order 
of S 2 ) from iv /(2b), but this difference can be neglected 
since it leads to the corrections of order of S 3 during the 
encoding/retrieval steps. 

By performing the straightforward calculations, one 
can arrive at the final answers for the performance pa- 
rameters st, sq, and s z : 



s = -85 2 , st = 1 — 4(5 , s z = 1 - 16<5 2 



(15) 



The deviation from ideal values is of order of 5 2 ~ l/N, 
which is about 5% for our simulations with TV = 20 (be- 
low, we will demonstrate the agreement in more detail). 

In our simulations, two types of initial conditions are 
considered: (i) the electronic spin is in the state = 1, 
Sx = Sy — 0, for evaluation of the s z performance mea- 
sure, and (ii) the electronic spin is in the state s. 



(i) 



0, for evaluation of the st and sq perfor- 



mance measures. Then we can use the values s z , st, 
and so to determine the fidelity of the protocol (see sec- 
tionlV|l . In order to simulate the initial state of the polar- 
ized bath, we take the bath's state as a random superpo- 
sition \r) of all possible basis states, apply the operator 
exp (—jl z ) to the state |r) (such an operator is easily 
implemented using the Chebyshev's polynomials expan- 
sion described in detail in Ref. ran , and normalize the 




FIG. 3: Final step of the protocol, retrieval. The values 
s x (a), st (b), so (c), and s z (d) are shown, obtained for 
AP — 0.2 with TV = 20. The solid curves represent the re- 
sults obtained with the bath's initial conditions obtained by 
application of the operator exp (— 7-P) to the random state, 
the dashed curves correspond to the analytical solution. 
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FIG. 4: Final step of the protocol, retrieval. The values s x (a), 
st (b), so (c), and s z (d) are shown, obtained for AP = with 
TV = 20. The solid curves represent the results obtained with 
the bath's initial conditions obtained by application of the 
operator exp (— 7-P) to the random state, the dashed curves 
correspond to the analytical solution. 



resulting state. Such a state corresponds to the bath's 
density matrix p — (l/Z) exp (— 7P). 

This construction of the bath's initial conditions can 
lead to statistical errors if the number of relevant states is 
not large enough. To ensure that the constructed initial 
condition is valid for TV = 20, we did the simulations 
with Ak — A, and compared them with the analytical 
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FIG. 5: The arrangement of the lattice of the bath spins 
and the electron density. The bath spins sites are marked 
by the solid circles, and the position of the maximum of the 
electron density is marked as <g>; the dashed line schematically 
shows the contour line where the electron density is equal to 
exp ( — 1/2) of its maximum value. 



results given in the previous section. The comparison is 
presented in Figs.EElfor AP = 0.2 (Fig.EJ) and AP = 
(Fig. 0J, for s z , so, st, and also for s x itself. Even in 
worst case of AP = 0.2, the simulation results remain 
close to the analytical ones, and at times close to the 
information retrieval time, the difference does not exceed 
~ 0.05. For smaller AP, the difference is smaller, and 
almost disappears for AP = 0, see Fig. 0] 

For the simulations, wc consider N = 20 bath spins 
which are placed at the sites of a 4 x 5 piece of a rect- 
angular lattice with the lattice constants a y and a x . We 
assume that the electrostatic potential inside the dot can 
be approximated as a 2-D harmonic well, so that the 
electron density in the dot is a 2-D Gaussian with the 
half- widths w x — (3/2)a x and w y — 2a y along the x- 
and y-axes. We set the maximum of the electron den- 
sity is slightly displaced with respect to the center of the 
lattice by l y — 0.2, l x =0.1 to prevent artificial symme- 
try effects. This arrangement is schematically shown on 
Fig. El where the bath spins' sites are marked by the solid 
circles, and the position of the maximum of the electron 
density is marked as ®; the dashed line schematically 
shows the contour line where the electron density is equal 
to exp (—1/2) of its maximum value. The values of Ak 
were taken proportional to the electron density with the 
proportionality factor 1, and were spread from 0.96 to 
0.31. 

The external field Hq in our simulations was taken 
equal to its optimal value, which compensates the Over- 
hauser field Ho v h = —PJ2k^-k/(2g e ^B) phis a small 
correction, H = -H 0vh — M 3 /(2M 2 g e ^B) (see Ea.lT2l. 
The memory protocol was simulated as follows: 

1. The bath's initial polarized state was obtained by 
acting with the operator exp (— 7P) to the random 
state; the initial state of the electron spin was | f). 
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FIG. 6: Performance of the memory protocol as a function 
of AP for w x = (3/2)a x and w y = 2a y . (a): the values of 
St (solid squares) and s z (open circles); (b): the values of so 
(solid circles). The values of st and s z are very close to each 
other, and the corresponding curves almost coincide on the 
graph. 
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FIG. 7: Performance of the memory protocol as a function of 
AP for w x — 3a x /(2\^2) and w y = \^2a y . (a): the values of 
st (solid squares) and s z (open circles); (b): the values of so 
(solid circles) . In spite of noticeably larger spread of the cou- 
plings At , retrieval values decrease only by 5-7% (depending 
on AP) in comparison with Fig. El 



2. The electron spin is allowed to interact with the 
bath, and the time when s z reaches its minimum is 
assumed to be the ejection time t e . 

3. The ejection of the first electron is described ac- 
cording to Eq. 00 

4. The second electron is injected in the state | !), 
so the compound system's density matrix has the 
form of Eq. 

5. The second electron is allowed to interact with the 
bath, and the time when s z reaches its maximum 
is assumed to be the retrieval time t r . 

The fidelity of the protocol is indirectly determined by 
the values of sr, s z , and sq at time t r : for ideal memory, 
sr(ir) and s z (t r ) would be one, and so(t r ) would be zero. 

The results for different values of AP are presented in 
Fig. Both s z and st are very close to each other, and 
both decrease approximately linearly with decreasing po- 
larization. Similarly, sq deviates from zero linearly with 
decreasing polarization. 

For comparison, in Fig. we present the results for 
the case of w x and w y decreased by a factor of 1/ \[2 with 
respect to Fig. (so that the dashed line of Fig. |S] corre- 
sponds now to the contour of the electron density equal 
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•so = -0.042, and s z = 1.0 - 0.084, and the numer- 
ical results are st — 1-0 — 0.028, so = —0.05, and 
s z = 1.0 — 0.082. The agreement between simulations 
and analytical results is reasonable, and we can consider 
our simulations to have a precision of 5-8%. 

V. MINIMAL FIDELITY OF THE MEMORY 
PROTOCOL 



FIG. 8: Performance of the memory protocol as a function of 
AP for the unoptimized external field, for w x = (3/2)a x and 
w y — 2a y . (a): the values of st (solid squares) and s z (open 
circles); (b): the values of so (solid circles). The maximum 
drop in performance in comparison with Fig. E] is only about 
6%, and is noticeable only for small AP. 



to e _1 of the maximum value). This corresponds to a no- 
ticeably larger spread of Ak, from 0.92 to 0.09. One can 
see that the perfromance of the protocol is determined 
by the value of s z . In spite of considerably larger spread 
of Ak, the change in performance measures s z ,sq,st is 
small, from 5%-7% depending on the bath polarization. 
This is expected: the smaller values of Ak give smaller 
contribution to the evolution, and, correspondingly, do 
not produce large errors in the course of the system's 
evolution. 

We note here that the quantum memory protocol de- 
scribed in the introduction implies fine adjustment of the 
external field. The optimal value of Hq includes the small 
correction M 3 /(2M 2 ) ~ l/N. While the Overhauser field 
can be measured rather precisely, the small correction 
Mz/{2M2) could be hard to determine in practice. How- 
ever, this is a l/N correction. According to our calcula- 
tions, omission of this correction does not affect much the 
performance of the scheme. Fig. [S] shows the results ob- 
tained for the same parameters as Fig. |H| except for the 
unoptimized external field (i.e., Hq — —Hovhi without 
the correction Ms/(2M2))- The maximum drop in per- 
formance is only 6% (for AP = 0), and goes to almost 
zero with increasing AP. 

Thus, our results show that the quantum memory pro- 
tocol remains feasible for realistic situations, where the 
polarization of the nuclear spins is incomplete, and the 
values of hyperfine couplings Ak are spread over a large 
interval. Also, the external field does not need to be ad- 
justed with precision of l/N ~ 10~ 4 for the quantum 
memory protocol to work. We believe that our results 
for N = 20 provide adequate estimates of the perfor- 
mance of the protocol for realistic dots with N = 10 4 . 
In order to understand the corrections associated with 
extending from N = 20 to larger N, let us compare the 
analytical estimates given in Eq. I|15|) with the numeri- 
cal results. For the case shown in Fig. El the analytics 
gives s T = 1.0-0.008, s = -0.016, and s z = 1.0-0.03, 
while the numerical simulations give st = 1.0 — 0.014, 
s = -0.022, and s z = 1.0 - 0.036. Similarly, for the 
case shown in Fig. [7] the analytics gives st = 1.0 — 0.021, 



We now show how the values of So, st and s z calculated 
above determine the fidelity of the memory protocol. The 
electron initially is in the state |0o) = a \ I) + P\ t)j 
so that the electron's reduced density matrix initially is 
T o = \4>o){<t>o\i an d the density matrix of the compound 
system is Rq — tq ® p, where p is the initial density 
matrix of the bath. At the end of the retrieval stage, the 
final state of the compound system is described by the 
density matrix Rt, so that the electron's reduced density 
matrix is Tf = Tr gRf, where Tr# is the trace over the 
bath spins. Ideally, we want tj — tq, but in reality the 
initial and the final density matrices of the electron are 
different. The fidelity for the protocol for the density 
matrix r is F Ta = Tr[T T/] = (0o|t/#o)- Accordingly, 
the minimal performance of the memory protocol (the 
fidelity, F) is defined in a standard way2£, as 

P = Mm 0o [^ o |r / |^o)], (16) 

which is minimized over all possible initial states \4>o). 

To calculate the fidelity, we use the ideas of the quan- 
tum process tomography^. First, let us analyze in de- 
tail the dynamics of the density matrix. At the encoding 
stage, the evolution of the compound system is described 
by the unitary operator U\ = exp [—iH.t e ], so by the end 
of the encoding stage, the state of the compound system 
is Pi = U\RqUI = U\{t (g) p)U\. During the ejection of 
the first electron and injection of another one, the evo- 
lution is not unitary (see Eas. I8I9|I : after injection of the 
second electron, the state of the compound system is: 

R e = P^P! + S-RiS + , (17) 

where P± = \ [){[ \ is the projection operator on the 
state | I) of the electron spin. At the retrieval stage, 
the evolution is again unitary, described by the operator 
U2 = exp [— i7i[t r — t e )] so the final state of the com- 
pound system is 

R f = UiP^iro^pplP^ 

+ U 2 S-U 1 {t g ® p)u\S + Ul (18) 

In order to analyze the evolution of the electron's re- 
duced density matrix, we consider the evolution super- 
operator C for the electron spin 

C[r }=r f = Tv B R f . (19) 

The superoperator C is linear. By representing the initial 
electron's density matrix in the form tq = (l/2)l+b x S x + 
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b y S y + b z S z , where 1 is the identity matrix 2x2, we 
obtain for the final density matrix Tf. 

T f = £{t } = r° f + b x T x + b y r) + b z r z f (20) 

where 

r° f = £[(1/2)1], rf=C{S% 

r) = C[S% rf=£[S z }. (21) 

Each of the matrices r9, . . . rf can be expanded in the 
basis 1, S x , S y , and S z , so we just need to determine 
the expansion coefficients. This problem is considerably 
simplified by the fact that the superoperator C is invari- 
ant with respect to the rotations of the compound sys- 
tem in the x-y plane, i.e. with respect to the operators 
exp [iaJ z ], where J z — S z + J^k ana - a IS ^ ne rotation 
angle. For example, let us apply the rotation by the an- 
gle 7r in the x-y plane (the operator exp[i7rJ z ]) to the 
matrix t°. By doing that, we obtain that 

Tr(rp x ) = -Tr(exp [-inJ z }T° f exp [iirJ z }S x ). (22) 

On the other hand, exp [— in J z )t® exp [iir J z ] ~ t9, so 
that Tt(t®S x ) — 0. By applying similar symmetry argu- 
ments, we derive the following results: 

T° f = (1/2)1 + S z u 

Tf = S Z W 

rf = S x v Q cos £ + S y v sin £, 
rj = S v vq cos £ — S x Vo sin £, 

where the last equation is obtained by applying 
the operator exp [z(7r/2) J z ] (rotation by ir/2 in the 
x-y plane) to the matrix rj, and noticing that 
exp[-i(7r/2)J z ]S' a: exp[i(7r/2)J z ] = S v . In Eqs. E2I u , 
wo, vq, and £ are numerical parameters characterizing 
the fidelity, which are related in a simple manner to the 
parameters sq, st, and s z obtained above from numerical 
simulations. We obtain s z by using the initial electron's 
density matrix tq = (1/2)1 + S z , so that s z = uq + Wq. 
Similarly, sq and st are obtained from the initial elec- 
tron's density matrix To = (1/2)1 + S x , so that sq = uq 
and st = vq. We note that the action of the superopera- 
tor C transforming S x and S y into and rj involves the 
rotation in the x-y plane by the angle £. As we explained 
above, this rotation can be, in principle, corrected, and 
we do not take it into account. Correspondingly, we ne- 
glect the reduction of fidelity caused by non-zero value 
of £, by putting £ = 0. 

By combining Eqs. |201 an d [23 the fidelity F To for a 
given initial density matrix To = (1/2)1 + b x S x + b y S v + 
b z S z can be expressed as 

F T0 = (1/2) [l + v + b z u + b 2 z {w - v )] , (23) 

where we used the fact that for initial pure state of the 
electron, b\ + b 2 = 1 — 6^, and we need to minimize the 



fidelity i12.'il) over all possible values of b z , i.e. over all b z £ 
[— 1 , 1] . The minimum value of this expression is achieved 
either at the extremum point b z o = — uq/[2(wq — vq)], or 
at the ends of the interval. Correspondingly, the minimal 
fidelity i*Mi n of the protocol is equal to the smallest of 
the three quantities 



_ 1 + s z - 2s 
F 2 - g ' 

^3 = ^(1 + st - s 2 /[A(s z -so- st)}) . (24) 

With this fidelity measure, our results described above 
and depicted in Figs. El an d El can be presented as 
a single graph, Fig. El One can see that the minimum 
fidelity is F Min - 75% for P = 80%, while as P -> 100%, 
-pMin —* 96%. The fact that the latter value is different 
from 100% is due to the modest number of bath spin N. 
Indeed, as seen from Eq. El for large N, the value of 
Fmiu for completely polarized bath differs from unity by 
the terms of order of 1/^V, and therefore goes to zero for 
large N. 

Using the analytical approach described in Sec. |ri] we 
can estimate the minimum fidelity for small AP in the 
approximation of all Ak = A = const. Indeed, the most 
significant contribution to s z , sq, and st is given by the 
states with / = — M (note that in our consideration 
M < 0); the contribution of the states with I = —M + n 
is of order (AP) n (where n is a positive integer). Thus, 
to calculate the values s z , sq, and st up to linear terms 
in AP, we take into account only two states, I = —M 
and / = — M + 1, which have relative statistical weight 
1 — AP/2 and AP/2, correspondingly, and perform av- 
eraging over all M according to Eq. El The calculations 
give the following results: 

A P 

s z = 1-— (5 + cos 4 7 o + sin 4 7 o) + 0((AP) 3 / 2 ) 
A P 

s ° = - f j-(5-cos 4 7 o-sin 4 7o ) + 0((AP) 3 / 2 ) 
A P 

st - 1-— (2-cos 70 ) + O((AP) 3 / 2 ) (25) 

where 70 = n/^/2, and from these values we obtain 

pyrin = 1 - 1.38386AP + 0((AP) 3/2 ) . (26) 

This line is shown in Fig. El as a dashed line. The dif- 
ference between this analytical result and the numerical 
data comes either from the assumption Ak = A = const 
used in the derivation above, and from the (AP) 3 / 2 cor- 
rections, which can be important for AP 0.1. 

Visible deviation of the fidelity from unity at 100% 
polarization results from finite size of the sample used 
in our simulations. Indeed, the analytical calculations 
of Ref. [Tl| and Section IIVI of this paper suggests that 
at 100% polarization the error arises from the inhomo- 
geneous hyperfine coupling, and hence scales as 1/ \/~N. 
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FIG. 9: Minimal fidelity Pnim of the memory protocol as a 
function of AP. Black circles: for w x = (3/2)a x and w y — 
2a y (same parameters as in Fig. Open circles: for w x — 
3a x /(2y/2) and w y = y/2a y (same parameters as in Fig. [7J. 
Open squares: unoptimized external field, for w x — (3/2)a x 
and w y = 2a y (same parameters as in Fig. 0. The dashed 
line corresponds to the analytical solution i'2(>ti for small AP. 



While it corresponds to few percent for the simulations 
presented in Fig.EI(iV = 20), for realistic dots (N = 10 6 ), 
it is well below 10~ 3 . 



fidelity of F ~ 80%. Fidelity is determined mostly by 
the retrieved value of s z , and increases linearly with in- 
creasing polarization of the bath. 

Since direct simulation of a real system with ~ 10 nu- 
clear spins is currently impossible, we consider the sys- 
tem of N = 20 spins coupled to the electron spin. Such a 
step requires justification, and accordingly a major part 
of our work is devoted to understanding the reliability 
of our simulations. By comparing different aspects of 
the numerical calculations with the analytical results, we 
show that our results give adequate estimates for realistic 
situations, with a precision of 5-8%. 

The errors considered in this work may be reduced by 
changes to the storage protocol or altered by consider- 
ing of more realistic initial conditions. For example, it is 
expected that at high polarizations, created by electron 
spin-nuclear spin transfer techniques (e.g., Ref. 0), the 
nuclear spins may be in states with favorable symmetry 
properties^ allowing for reduced error due to polariza- 
tion. Similarly, estimating the current value of the nu- 
clear spin polarization reduces the error associated with 
the uncertainty in polarizationiS*^. More detailed work 
will be necessary to consider the effects of these complex- 
ities. 



VI. CONCLUSIONS 

We study numerically the long-lived memory for elec- 
tronic spin in a quantum dot taking into account the 
errors inherent in the protocol for a realistic system: the 
incomplete polarization of the bath of nuclear spins and 
the spread in the hyperfine interactions between the elec- 
tron and the nuclei in the quantum dot. Although the 
imperfections deteriorate the fidelity of the quantum in- 
formation retrieval, the scheme remains feasible for the 
spin bath polarizations exceeding 80%, with a minimum 
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